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Abstract 



We present updated predictions for the total cross section for Higgs boson production 
by gluon-gluon fusion in hadron collisions. Our calculation includes the most advanced 
theoretical information available at present for this observable: soft-gluon resummation 
up to next-to-next-to-leading logarithmic accuracy, the exact treatment of the bottom- 
quark contribution up to next-to-leading order, and two-loop electroweak effects. We 
adopt the most recent parametrization of parton distribution functions at next-to- 
next-to-leading order, and we evaluate the corresponding uncertainties. In comparison 
with our previous central predictions, at the Tevatron the difference ranges from +9% 
for mjj = 115 GeV to —9% for nin = 200 GeV. At the LHC the cross section is instead 
significantly increased. The effect goes from +30% for m# = 115 GeV to +9% for 
m H = 300 GeV, and is mostly due to the new parton distribution functions. We also 
provide new predictions for the LHC at y/s = 10 TeV. 
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The Higgs boson is a key ingredient of the Standard Model (SM), but it has so far eluded 
experimental discovery. Direct searches at LEP lead to a 95 % CL lower limit of win > 114.4 GeV 
on the mass m# of the SM Higgs boson pQ . At the LHC the Higgs boson can be discovered over 
the full mass range up to mn ~ 1 TeV within a few years of running. At the Tevatron, the CDF 
and DO experiments are now becoming sensitive to a Higgs signal at mn ~ 170 GeV [2]. 

The dominant mechanism for SM Higgs boson production at hadron colliders is gluon-gluon 
fusion, through a heavy-quark (mainly, top-quark) loop. When combined with the decay channels 
H — > 77 and H — > ZZ , this production mechanism is one of the most important for Higgs boson 
searches and studies over the entire range, 100 GeV ^itlh ^ 1 TeV, of Higgs boson mass to be 
investigated at the LHC. In the mass range 140 GeV ^ 180 GeV, gluon fusion, followed by 
the decay H — > WW — > lulu, offers the main discovery channel of the Higgs boson at the LHC 
and also at the Tevatron, thanks to the strong angular correlations of the charged leptons in the 
final state. 

In QCD perturbation theory the leading order (LO) contribution to the gg — > H cross section 
is proportional to a|, as being the QCD coupling. The QCD corrections have been computed 
at next-to-leading order (NLO) [21 H] in the heavy-top limit, and with full dependence on the 
masses of the top and bottom quarks (4]. Next-to-next-to-leading order (NNLO) corrections 
have been obtained in the heavy-top limit [5]. These QCD corrections, which are dominated by 
radiation of soft and virtual gluons [6], lead to a substantial increase of the LO result. The QCD 
computation up to NNLO has been consistently improved by adding the resummation of soft- 
gluon logarithmic contributions, up to next-to- next-to-leading logarithmic (NNLL) accuracy [TJ. 
The ensuing NNLL+NNLO results are nicely confirmed by the more recent computation [8]-[TU] 
of some of the soft-gluon terms at N 3 LC0. 

In this Letter we present an update to the NNLL+NNLO results of Ref. [TJ. We refrain from 
repeating the theoretical details of the calculation, which can be found in Ref. [TJ, and we only 
describe the main changes implemented in the present work. In our previous analysis the bench- 
mark predictions were obtained using the MRST2002 NNLO partons [TT] . Since then, there have 
been important modifications in the extraction of parton distribution functions (PDFs). The cal- 
culation of the NNLO splitting functions^ was completed |12j . Furthermore, a more sophisticated 
treatment of heavy-quark thresholds has been introduced in [T3] , resulting in a significant modi- 
fication of the gluon and light quark densities in the recent MSTW2008 set [15]. For example, at 
x ~ 0.01, relevant for a Higgs boson of ~ 120 GeV at the LHC, the gluon distribution increases 
by about 6% with respect to the MRST2002 fit. The value of as(mz) also had a non-negligible 
change from 0.1154 to 0.1171. Considering that the total cross section is completely dominated 
by the gluon-gluon fusion channel, and that the lowest order contribution starts at 0(a^), with 
sizeable corrections at higher orders, it is not surprising that the mere change from MRST2002 to 
MSTW2008 partons can result in an increase of more than 10% in the production cross section, 
making an update mandatory. The change in the PDFs does not result in such a dramatic increase 
of the cross section at the Tevatron, since at x ~ 0.06, relevant for the production of a Higgs boson 
of mu ~ 120 GeV, the gluon distribution is reduced by about 4%, but the decrease is partially 
compensated by the increase of the partonic cross section due to the larger coupling constant. 

^Note, however, that the perturbative information available at present does not allow to consistently extend the 
computation to N 3 LL accuracy. 

■I-Thc MRST2002 fit was based on approximate expressions for the evolution kernels [13]. 
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Besides the important effect of the PDFs, there are other theoretical reasons for revisiting 
the computation of the Higgs cross section at hadron colliders. In particular there has been 
an important effort to evaluate the electroweak (EW) corrections arising from W and Z boson 
coupling to the Higgs and to both light and heavy quarks in the loop [16]. The recent computation 
of Ref. p2] takes into account those contributions by avoiding the complications in the two-particle 
threshold using the complex-mass scheme. The EW corrections turn out to be of the order of a 
few percent, with a sign depending on the Higgs mass. The main uncertainty in the EW analysis 
comes from the fact that it is not completely clear how to take them into account in practice. In 
the partial factorization scheme of Ref. [TTj the EW correction applies only to the LO result. In 
the complete factorization scheme instead, the EW correction multiplies the full QCD corrected 
cross section. Since QCD corrections are sizeable, the latter choice has a non-negligible effect 
on the actual impact of EW corrections in the computation. The recent analysis of higher-order 
QCD and EW corrections presented in Ref. |18j . performed on the basis of an effective Lagrangian 
approach, supports the complete factorization hypothesis, suggesting that EW corrections become, 
to a good approximation, a multiplicative factor of the full QCD expansion. 

The predictions we present below are obtained as follows. We first consider the top-quark 
contribution in the loop, and perform the calculation up NNLL+NNLO in the large-mi limit. 
The result is rescaled by the exact m t dependent Born cross section, since this is known to be an 
excellent approximation for the top-quark contribution. This resummed top-quark contribution 
provides the bulk of the Higgs cross section at hadron colliders. We then consider the bottom- 
quark contribution (more precisely, the bottom contribution and the top-bottom interference). 
Since in this case the effective theory approach is not applicable, we follow Ref. [IS] and we 
include this contribution up to NLO only (but still computed with NNLO MSTW2008 partons), 
by using the program HIGLU [I]. Finally, we correct the result by including the EW effects 
evaluated in Ref. [17] in the complete factorization hypothesis. We set the heavy-quark masses 
to m t = 170.9 GeV and m b = 4.75 GeV, the latter consistently with the MSTW2008 set. Our 
central predictions (<r best ) are obtained by setting the factorization ([if) and renormalization (fin) 
scales equal to the Higgs boson mass. 

Our results for the Tevatron at y/s = 1.96 TeV and the LHC at \/s = 10 TeV and \/s = 14 
TeV are presented in Tables 1, 2 and 3, respectively. Comparing to our previous predictions (see 
Tables 1 and 2 of Ref. [7]), the cross sections change significantly. At the Tevatron the effect 
ranges from +9% for ran = 115 GeV to —9% for ran = 200 GeV. At the LHC the effect goes from 
+30% for rrtu = 115 GeV to +9% for ran = 300 GeV. It is worth noticing that at the LHC more 
than half of the increase arises from the modification in the gluon distribution and the coupling 
constant. 

The bottom contribution, dominated by bottom-top interference, is small and negative. The 
different treatment of this contribution with respect to the previous analysis [7] results in an 
increase of the cross section from about 7% (m# = 115 GeV) to 4% (m# = 200 GeV) at the 
Tevatron and from 9% (m H = HO GeV) to 2% (m H = 300 GeV) at the LHC. The inclusion of the 
EW corrections results in an increase of the cross section by about 5% for mu ^ 160 GeV, and a 
decrease by about 2% for 200 GeV £ m H £ 300 GeV. 

Our results for the Tevatron can be compared to those recently presented in Ref. [18], obtained 
using the same set of PDFs. Besides the different choice for the bottom-quark mass and the 
implementation of an effective Lagrangian calculation for the EW corrections, the main difference 
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with our work arises in the calculation of the top-quark contribution to the cross section. In 
Ref. [18] the latter contribution is computed up to NNLO but choosing p, F = fi R = m H /2, as 
an attempt to mimic the effects of soft-gluon resummation beyond NNLO. The final numerical 
differences at the Tevatron turn out to be small and of the order of a few per mille at the lowest 
masses, increasing to 2.5% at m# = 200 GeV. 

We now discuss the various sources of uncertainty affecting the cross sections presented in 
Tables 1, 2 and 3. The uncertainty basically has two origins: the one coming from the partonic 
cross sections, and the one arising from our limited knowledge of the PDFs. 

Uncalculated higher-order QCD radiative corrections are the most important source of uncer- 
tainty on the partonic cross section. A method, which is customarily used in perturbative QCD 
calculations, to estimate their size is to vary the renormalization and factorization scales around 
the hard scale m#. In general, this procedure can only give a lower limit on the true uncertainty. 
The uncertainty is quantified here as in Ref. [7]: we vary independently fip and in the range 
0.5m# < Hf,Hr < 2m#, with the constraint 0.5 < Hf/hr < 2. The results are reported in Tables 
1,2 and 3. The scale uncertainty is about ±9 — 10% at the Tevatron and ranges from about ±10% 
(m H = 110 GeV) to about ±7% (m H = 300 GeV) at the LHC (y/s = 10 and 14 TeV). These re- 
sults are consistent with those of Ref. [7J; in particular, we note that the effect of scale variations 
in our resummed calculation is considerably reduced with respect to the corresponding NNLO 
result. At NNLO the scale uncertainty is about ±14% at the Tevatron and ranges from about 
±12% (m H = HO GeV) to about ±9% (m H = 300 GeV) at the LHC (y/s = 10 and 14 TeV). The 
reduction is more sizeable at the Tevatron, where the resummation effect is more important. 

Another source of perturbative uncertainty on the partonic cross sections comes from the 
implementation of the EW corrections. Our results are obtained in the complete factorization 
scheme discussed above. The partial factorization scheme would lead to a change of our results 
ranging from about —3% (m# = 115 GeV) to ±2% (mjj = 200 GeV) at the Tevatron and from 
about -3% (m H = 110 GeV) to ±1% (m H = 300 GeV) at the LHC. 

A different source of perturbative uncertainty comes from the use of the large-m^ approximation 
in the computation of the partonic cross section beyond LO. The comparison between the exact 
NLO cross section and the one obtained in the large-m t approximation (but rescaled with the 
full Born result, including its exact dependence on m t ) shows that the approximation works well 
also for m>H > m t . This is not accidental: the higher-order contributions to the cross section are 
dominated by relatively soft radiation, which is weakly sensitive to the mass of the heavy quark in 
the loop at Born level. This feature persists at NNLO and thus it is natural to assume that having 
normalized our resummed result with the exact mi-dependent cross section, the uncertainty due 
to the large-mt approximation should be of the order of few percent, as it is at NLO. The effect 
of finite- m t corrections is discussed in Refs. [191 120] . 

The other important source of uncertainty in the cross section is the one coming from PDFs. 
Modern PDF sets let the user estimate the experimental uncertainty originating from the accuracy 
of the data points used to perform the fit. The most recent MSTW2008 NNLO set provides 40 
different grids that allow us to evaluate the experimental uncertainties according to the procedure 
discussed in Ref. [UJ. The outcoming uncertainties (at 90% CL) are reported in Tables 1,2 and 
3. At the Tevatron the effect ranges from ±6% (m# = 115 GeV) to about ±10% (m# = 200 
GeV), while at the LHC it is about ±3% (±3 — 4% at yfs = 10 TeV) in the mass range we have 
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considered. We note that at the LHC the PDF uncertainty is quite small, and, in particular, it is 
substantially smaller than the uncertainty from missing higher-order perturbative contributions, 
as estimated from scale variations. On the contrary, the PDF uncertainty at the Tevatron is larger. 
This is a consequence of the fact that the Tevatron probes relatively larger values of x, where the 
gluon density is less constrained. By using the MRST2006 set [Hj our result for m# = 170 GeV at 
the Tevatron would be a best ( 2006 ) = 0.395±g;^f (scale) (PDF), which is 13% smaller than the 
one reported in Table 1, cr best = 0.3491qq27 (scale)!^'^ (PDF), and marginally compatible with 
it. 

We finally point out that a related and important uncertainty is the one coming from the value 
of the QCD coupling. In modern PDF sets ces( m z) is obtained together with the parton densities 
through a global fit to the available data, and thus there will be a correlation between the error on 
a s{ m z) and that on the gluon density. Since the gluon fusion process starts at O(o.g), it is easy 
to foresee that the uncertainty on a>s(mz) may have an important impact on the cross section. 
Neglecting correlations with the gluon density, a 3% uncertainty on a s (m^) would lead to an 
effect of about ±9 — 10% on the production cross section at both the Tevatron and the LHC. 

To summarize, we have presented updated predictions for the cross section for Higgs boson 
production at the Tevatron and the LHC. The results are based on the most advanced theoretical 
information available at present for this observable, including soft-gluon resummation up to NNLL 
accuracy and two-loop EW corrections. In comparison with the central predictions of Ref. [7J, 
at the Tevatron the difference ranges from +9% to —9% for 115 GeV ^m#^ 200 GeV. At the 
LHC the effect goes from +30% to +9% for 115 GeV ^m H ^ 300 GeV, and is mostly due to the 
new PDFs. We have then reviewed [7J the uncertainties that affect the Higgs production cross 
section, and we have shown that they are still relatively large, especially at the Tevatron. The 
above uncertainties should be taken into account in Higgs boson searches and studies at both the 
Tevatron and the LHC. 

We thank Babis Anastasiou and Christian Sturm for useful discussions. We are grateful to 
Stefano Catani for helpful discussions and comments on the manuscript. The work of D.deF. 
was partially supported by ANPCYT, UBA-CyT and CONICET. D.deF. is grateful to the US 
Department of Energy (Contract No. DE-AC02-98CH 10886) for providing the facilities essential 
for the completion of his work. 
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m H 


^-best 


Scale 


PDF 


100 


1.861 


+0.192 
-0.174 


+0.094 
-0.101 


105 


1.618 


+0.165 
-0.149 


+0.085 
-0.091 


110 


1.413 


+0.142 
-0.127 


+0.077 
-0.083 


115 


1.240 


+0.123 
-0.110 


+0.070 
-0.075 


120 


1.093 


+0.107 
-0.095 


+0.065 
-0.069 


125 


0.967 


+0.094 
-0.083 


+0.059 
-0.063 


130 


0.858 


+0.082 
-0.072 


+0.054 
-0.058 


135 


0.764 


+0.073 
-0.063 


+0.050 
-0.053 


140 


0.682 


+0.065 
-0.056 


+0.046 
-0.049 


145 


0.611 


+0.057 
-0.049 


+0.042 
-0.045 


150 


0.548 


+0.051 
-0.044 


+0.039 
-0.042 


155 


0.492 


+0.045 
-0.039 


+0.036 
-0.038 


160 


0.439 


+0.040 
-0.034 


+0.033 
-0.035 


165 


0.389 


+0.036 
-0.030 


+0.030 
-0.032 


170 


0.349 


+0.032 
-0.027 


+0.028 
-0.029 


175 


0.314 


+0.029 
-0.024 


+0.026 
-0.027 


180 


0.283 


+0.026 
-0.021 


+0.024 
-0.025 


185 


0.255 


+0.023 
-0.019 


+0.022 
-0.023 


190 


0.231 


+0.021 
-0.017 


+0.020 
-0.021 


195 


0.210 


+0.019 
-0.015 


+0.019 
-0.020 


200 


0.192 


+0.017 
-0.014 


+0.018 
-0.019 



Table 1: Cross sections (in pb) at the Tevatron (fip = [ir = mu) with y/s = 1.96 TeV, 
using the MSTW2008 parton densities. 
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^best 


Scale 


PDF 


m H 


0-best 


Scale 


PDF 


m H 


0-best 


Scale 


PDF 


100 
110 
120 
130 
140 
150 
160 


44.12 
36.99 
31.48 
27.11 
23.58 
20.69 
18.07 


+4.24 
-4.44 

+3.43 
-3.60 

+2.83 
-2.96 

+2.35 
-2.48 

i 1 no 

-2.10 

+1.69 
-1.80 

+1.44 
-1.53 


+1.07 
-1.39 

+0.88 
-1.14 

+0.75 
-0.96 

+0.64 
-0.82 

i n rift 
-|-U.OO 

-0.71 

+0.50 
-0.62 

+0.44 
-0.55 


170 
180 
190 
200 
210 
220 
230 


15.63 

13.78 

12.20 

10.97 

9.98 

9.14 

8.42 


+1.22 
-1.30 

+1.05 
-1.12 

+0.91 
-0.97 

+0.80 
-0.86 

-|-U. ' ^ 

-0.77 

+0.64 
-0.69 

+0.58 
-0.63 


+0.39 
-0.48 

+0.35 
-0.42 

+0.32 
-0.38 

+0.29 
-0.35 

i n o7 
-0.32 

+0.26 
-0.30 

+0.24 
-0.28 


240 
250 
260 
270 
280 
290 
300 


7.81 
7.29 
6.83 
6.44 
6.11 
5.83 
5.61 


+0.53 
-0.58 

+0.49 
-0.53 

+0.45 
-0.49 

+0.42 
-0.46 

i n /in 
-0.43 

+0.37 
-0.40 

+0.37 
-0.38 


+0.23 
-0.26 

+0.22 
-0.25 

+0.21 
-0.24 

+0.21 
-0.23 

i n on 
-0.22 

+0.20 
-0.22 

+0.19 
-0.21 



Table 2: Cross sections (in pb) at the LHC (\xp = = mu) with ^/s = 10 TeV using 
the MSTW2008 JIf parton densities. 
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^.best 


Scale 


PDF 


m H 


^.best 


Scale 


PDF 


m H 




Scale 


PDF 


100 
110 
120 
130 
140 
150 
160 


74.58 
63.29 
54.48 
47.44 
41.70 
36.95 
32.59 


+7.18 
-7.54 

+5.87 
-6.20 

+4.88 
-5.18 

+4.12 
-4.38 

+3.47 
-3.75 

+3.02 
-3.24 

+2.60 
-2.79 


+1.86 
-2.45 

+1.54 
-2.02 

+1.30 
-1.70 

+1.12 
-1.45 

+0.97 
-1.25 

+0.85 
-1.10 

+0.73 
-0.97 


170 
180 
190 
200 
210 
220 
230 


28.46 
25.32 
22.63 
20.52 
18.82 
17.38 
16.15 


+2.22 
-2.39 

+1.92 
-2.08 

+1.68 
-1.83 

+1.49 
-1.63 

+1.34 
-1.47 

+1.22 
-1.33 

+1.11 
-1.22 


+0.65 
-0.84 

+0.58 
-0.74 

+0.52 
-0.66 

+0.48 
-0.60 

+0.45 
-0.55 

+0.42 
-0.51 

+0.39 
-0.48 


240 
250 
260 
270 
280 
290 
300 


15.10 
14.19 
13.41 
12.74 
12.17 
11.71 
11.34 


+1.03 
-1.12 

+0.95 
-1.04 

+0.88 
-0.97 

+0.83 
-0.91 

+0.78 
-0.86 

+0.74 
-0.82 

+0.71 
-0.78 


+0.37 
-0.45 

+0.36 
-0.43 

+0.35 
-0.41 

+0.33 
-0.39 

+0.33 
-0.38 

+0.32 
-0.37 

+0.32 
-0.36 



Table 3: Cross sections (in pb) at the LHC (\xp = = mu) with ^/s = 14 TeV using 
the MSTW2008 parton densities. 
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